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ABSTRACT 

The primary problem in the calculation of the differential 
scattering cross section for hydrogen molecular ions is the derivation 
of an approximate expression for the probability that the hydrogen 
molecular ion will be oriented, prior to impact, in such a way as to 
avoid disruption during scattering. Using a simple model, an approxi- 
mate expression for the orientation probability is obtained. 
Differential cross sections for the scattering of hydrogen molecular 
ions from carbon and krypton are calculated using classical scattering 
theory and the derived orientation probability. Impact energies of 
100 to 1000 electron volts are considered. In addition differential 
cross sections for protons scattered from carbon and krypton are 


presented. 
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1. INTRODUCTION 

Bates, Cook, and Smith, have presented a classical theory of 
ion-molecular rearrangement collisions at high impact ee ot 
Using proton-methane collisions as an example a brief outline of the 
theory follows. 

In Figure 1 m, is the impacting proton, m: is one of the 
hydrogen atoms of the methane molecule and M. is the carbon atom of 
the methane molecule. In this collision m, =M < M- For capture 
to occur Vw We . From these two conditions and the conservation of 
of linear momentum and energy it can be shown that @5 ww 90° and 
6; w5 . Thus, for this case, capture of the hydrogen atom of the 
methane by the incident proton to form hydrogen molecular ions may 
occur at an angle of approximately 45°. Bates, Cook and Smith have 
calculated a cross section for the production of hydrogen molecular 
ions by this aroseesttl. 

An experimental difficulty may arise in a measurement of the 
cross section for the production of He by this double scattering 
process if the incident proton beam contains some significant number 
of Hp . If the cross section for scattering of Hp from the carbon 
atom is much larger tha the cross section for H, production by capture, 
then a measurement of the production cross section would be impos- 
sible unless a pure source of protons were used. In order to 
determine the relative magnitudes of these two differential scatter- 


ing cross sections a calculation of the differential cross section 


of i ions scattered from heavy atoms is presented. 
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2. THE MODEL 

The hydrogen molecular ion is linear with an inter atomic spacing 
of approximately two angstroms. For randomly oriented impacting HE 
ions this implies a possible difference in impact parameter for the 
individual atoms of as much as two angstroms. If the difference in 
impact parameters is great enough, the velocities of the two atoms 
after scattering may differ by enough to cause their relative energy 
to exceed the binding energy of the molecule. The molecule will 
then be disrupted and can not be considered an He scattering event. 

Even at zero difference in impact parameter disruption of the 
molecule may occur because of the following consideration. Ina 
spherically symmetric potential, such as that used in this calcu- 
lation, each particle remains in a plane, defined by the scattering 
axis and a line from the scatterimg axis to the particle, during 
and following the collision. If the angle between the planes is 
great enough the difference in the velocities of the two atoms, 
after scattering, may be such that the molecule is disrupted. 

Figure 2 shows the relationship between the scattering planes 
and the scattering axis. C@is the angle between the scattering 
planes. 

The condition that the molecule remains intact may be stated as 
D> x (Cv) where D is the binding energy, CV is the relative 
velocity of the two particles and 4% u(CVv)* gives the relative energy 
of the two particles. Using the symbol AV for the maximum allowable 
CV the Limiting condition for nondisruption becomes: 


D = } uw (Av). 
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Thus both the difference in impact parameter €p and the angle between 
the scattering planes must be such that the above condition is 
satisfied. This limitation on the orientation of the impacting 
molecule leads to an orientation probability that must be considered 
in the calculation of the differential scattering cross section for 
hydrogen molecular ions. 

For high impact energies and small impact parameters each atom 
of the He may be treated as a proton. That is, the effective charge 
may be taken as one for each particle. The calculation of the differ- 
ential cross section for He - ions may then be approached in the follow- 
ing way. For a given energy and impact parameter calculate the 
scattering angle, and the differential cross section for that angle, 
for one proton whose velocity is V,, the velocity of the H, . Then 
determine the probability that the "second" particle is oriented with 
respect to the "first" particle in such a way that its scattering 
angle is within A @ of 9 and its scattering plane makes an angle of 
less than A @ with the scattering plane of the first particle. The 
differential cross-section for He scattering will then be given by: 
dys (8) = 0,,+(6) P, where P is the orientation probability. This is 
then done over a range of impact parameters and for several energies. 
3, CALCULATION of 9,,(6) 

The calculation of o,,(6) for the "first" particle was done using 
the method suggested by Bates, Cook and smith’ 4 | The potential used 
for the calculation was a shielded potential of the forn, 

V(R) = - (Men New"+ ne a) due to pyatte 2 | Calculations were made 
for both carbon and krypton. Byatt lists several potentials in 


reference 2, all of which may be conveniently used in the computer 
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program in appendix 2. The cross section was calculated for impact 
parameters from 0 - 5 a, in increments of 0.05a, (a, is the radius of 
the Bohr orbit) and for energies from 100ev to 1000ev in steps of 
100ev. .The details of the calculation are shown inappendix l. 

4.. CALCULATION OF THE ORIENTATION PROBABILITY. 


The condition AV “(i )’ may, for convenience, be put in terms 
My 


of'the energy of the Hp. Bt = 45 (2M) Vp = M,V5 or My = te then 


AV = [4p a This limiting condition implies a relationship 
ce 


between €9@ and €® such that AV= F(C@, C%). CO is related to Co by 

a = oCe/oy0) SIN@ in the center of mass coordinate system. 

Therefore AV = G(C®,Co) gives a relationship between €@ and €op. Since 
Aree a constant for a given energy and impact parameter, there exists 
rome relationship between Co and €@ such that C@= £C€p). CO= £ (Co) 
derines a curve in a plane perpendicular to the scattering axis. The 
second particle must pass through the area on the plane bounded by 
this curve if disruption is to be avoided. 

Te arti y oriented He can be considered as defining a Senee 
the center of which is the "first" atom and whose radius is £, the 
inter atomic spacing. , 

“Th Figure 3, A is the area on the surface of the hemisphere con- 
tained within the surface defined by Co = f£(€@). The ratio of this 
ae to the area of a hemisphere of radius 2 gives the fraction of 
the He ions which are oriented in such a way as to avoid disruptions. 
Thus the orientation probability, P is given by: P = A/2mg*. No exact 
Poiteten has been obtained for the orientation probability. 


As an approximation to the orientation probability the following 


simplifications were made. 
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Assuming the target atom to be infinitely massive, there will be 
no chamge in the magnitude of the velocity of either atom of the 
hydrogen molecular ion. Since they then have the same initial speed 


V, When C® = 0 the particles are scattered in the same 


QO 3 


MW = Va = V,. 


plane as shown in Figure 4a. Then, assuming 49 is small 


a 
AV = VoL49 or 4G “ie) . When €@ = 0 (which implies Cp = 0) 
2 


the relationship between the individual atoms is as shown in Figure 4b. 
ts L 
Then AV = (V, SING) AD= 4D y or 4b = 4D \? _1_\ these approximations 
ee a Sit: SINO 
then give independent maxima for €® and CQ. AQ is related to Do by /. 


Lo 7 9) SING £9 or Lo a +6) SIN@ AD 4 
p p i 
“HG 
The above restriction on Ao implies that the second patticle must 


be contained between two coaxial cylinders whose axis is the scattering 


axis and whose radii are 9 - 4o and 9p + 4o respectively. The condition 


LD w 4D \% ob ae 
CHE SING 


puts a limit on the angular separation of the scattering planes. The 
limits on AD and Ap then confine the second particle to lie on the 
surface of a hemisphere of radius £ contained between the two scattering 
planes and the two cylindrical surfaces. Figure 5 illustrates. 

This area, on the surface of the hemisphere, divided by the area 
of the hemisphere then gives an approximation to the orientation 
probability. 

As a further approximation, this area was assumed small and taken 
to be the plane area given by Aw p(2A%) (240); which then gives, as 


an eatimate of the orientation probability, Pw 2040 CC 
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5. RESULTS AND BISCUSSION 
Everhart, et al, have calculated the differential cross section 
L3] 


for scattering from a screened Coulomb potential - The results are 


presented in terms of a, the screening length, given by: 
ae ae Ba 


and b, the collision diameter, given by: 
2 


where pp = the reduced mass of the system, 


<j 
tl 


the relative velocity of projectile and target. 

As a check on the calculation of the cross section for proton 
scattering the computer program used for the calculation was run using 
the screened Coulomb potential of carbon and a proton impact energy 
such that the ration b/a was equal to 1. The results of this calcula- 
tion were within 2.34% of those of Everhart et al. 

The actual calculations were made using the Byatt potentials of 
carbon and krypton with proton impact energys from 50 to 500 electron 
volts. The results are presented in the graphs of Figures 6 and 7. 

The cross section for the scattering of the hydrogen molecular 
ion for impact energies of 100 to 1000 electrons volts (corresponding 
to proton energies of 50 to 500 electron volts) then is the cross 
section for proton scattering multiplied by the orientation probability. 
These cross sections are presented in the graphs of Figures 8 and @. 

Under the assumption that the speeds of the two particles are the 
same following scattering, the AQ used in the orientation probability 


approximation is due only to the difference in the directions of the 
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0 10 120 30 0 D0 60 0 80 90 
Scattering angle-degrees (lab. coords.) 
Figure 6. Graph of the logarithm of the differential scattering cross 
section versus scattering angle for protons on carbon. 
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Figure 7. Graph of the logarithm of the differential scattering cross 
section versus scattering angle for protons on krypton. 
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0 10 20 30 40 50 60 70 80 90 
Scattering angle-degrees (lab. coords.) 
Figure 8. Graph of the logarithm of the differential scattering cross 
section versus scattering angle for hydrogen molecular ions on carbon. 
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Figure 9. 
section versus scattering angle for hydrogen molecular ions on krypton. 
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Graph of the logarithm of the differential scattering cross 
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velocity vectors. The speed of a particle scattered through an angle 


9 is given by: 


. 
V=V, i <ofntan (1- a) 


( +M 


where V = the speed of a projectile particle scattered through 
an angle § in the C.M. Coordinate systen, 
a Mass of the projectile particle, 
M.= Mass of the target particle 
and V,% initial speed of the projectile particle. 
Lf the speeds of particles scattered through angles 9 and § + C@ 
differ by enough that €V,, resulting from this difference, is greater 
than AV, the maximum for non disruption, then the Hy will be disrupted 
due to the speed difference alone. If 10,>>09, then AQ (ue to the 
difference in direction only) will be the more stringent condition and 
the error introduced by ignoring the difference in speeds may be assumed 
small. 
2M M | : : -. 
Thestacter  t » for scattering from carbon, is approximately 
mea No 
(4, ) 


meee . The factor (l1- cosQ§) is always less than one for scattering 
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angles between zero and 90 degrees. Thus the expression for V above 


can be expanded in a series. Takimg only the first two terms gives, 


2 (1- cos@) 


V= V, See 
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iL 
Therefore AV = /4D \ * , which is just the 49 due to the differ- 
V. - 
ence in direction Of the velocity vectors. 


+M,\ * 
Therefore 49 = t . For carbon, then eos 214 AW. 


MM sing 
This implies that the He will be disrupted due to the difference in the 
directions of the velocities of the two particles at a A9 which is much 
less than the 49, required to cause a disruption due to a speed differ- 
ence only. 

Because there is some difference in the speeds of the two particles 
the relative energy of the particles is somewhat greater than that 
predicted with the speeds assumed equal. Thus disruption will occur 
more often than predicted by the model used for the calculation. 
Therefore, the assumption of equal speeds tends to make the probabil- 
ity factor somewhat too large. 

The use of separate maxima for 46 and A® introduces an error in 
that the final probability expression contains the product of the 
maximum values, each of which was derived under the assumption that 
the other was zero. This error tends to make the probability factor 
too large. 

In view of the small values of P, (representative values are 
listed in Table 1.) the assumption that the area on the surface of the 
hemisphere can be approximated by a plane area would seem to be 
reasonable for, at least, the higher angles and energies. 

The conclusion is that the probability factor is in general some- 


what too large thus causing the calculated cross sections for 
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Energy 
(in ev) 


100 
100 
100 


500 


Table l. 


@ (Lab) 


32,0) 
45.62 
10.68 
88.88 
42.35 
10.37 
735.20 
47.62 


INO RSS 15) 


Carbon 
P 


1.95xl0 
Pata) 
2.10x10— 
9.18x10— 
7.08x10— 
2.30x10— 
1.72x10— 
2 107 


7.46x10- 
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Krypton 
Q (Lab .) P 
87.88 5.08x10— 
Mg 25 2 22a 
10 32 6.38x10 
83.45 2.68x10 
45.39 1.66x10— 
10.10 7.00x10- 
86.47 6.04x10_ 
46.80 3a] 2s 
10.13 2.32x10- 


Representative values of the orientation probability. 
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hydrogen molecular ion scattering to be too large. 

In Figure 10 the curve marked Oys(45° ) is the differential scatter- 
ing cross section for He from carbon at a scattering angle of 45° as 
a result of this calculation. The impact energies plotted are from 
100ev to 1000ev. The curve labeled o (He+CH, —» CHj3+H,) is the cross 
section for the capture of a hydrogen atom from the methane molecule 
by an impacting proton. This curve is a result of calculations based 
on the reaction cross section as derived by Bates, Cook and smitht 1 | 
Figure 10 shows the cross section for scattering of Hp from carbon at 
45 to be somewhat more than an order of magnitude greater than the 
cross section for the reaction H’ +CH, —» CHa +Hp . As a result, in 
order to measure the reaction cross section it would appear that a 


proton beam free of He will be necessary. 
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Figure 10. 


Comparison of reaction and scattering cross 
sections. 
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APPENDIX 1 


CALCULATION OF (j,+(9) 


The classical differential scattering cross section in the center 


of mass coordinate system is given by: 


Oe) = @ [sielggl]™ 1) 


where @is the impact parameter and@is the scattering angle in the 


center of mass coordinate system. 6is related to@by: 6=1M-2CU (2) 


i>) 


where U r. 
= ray 
ryi- $2, -- |" (3) 


in whichwis the largest positive root of 
2 
|- fan - FE (4) 
and A{ = the reduced mass of the system 
Vir) = the scattering potential. 
V. = the velocity of the impacting particle in the lab 


coordinate system when the target particles are considered 


EO be at rest. 


The potential as given by Byattt24 is 
2 r 
V(r) = Za (De“+Eel+ re (5) 
after converting V(r) to electron volts and letting 2fq, there results: 
An oc R PR oR 6 
V(R) = & (Ae + BOM + ce"*) (6) 
Where For Carbon For Krypton 
A= 204.136 328059 
B= -71.856 588.600 
= 31.029 63.7bo 
O= =12700 -1.086 
= -2.895 -4.961 
y= -8.808 ~26.110 


The quantity av, can be put in terms of the energy of the H,” 


in electron volts in the lab. coordinate system. 





2 = ]. Ext - oft + E yt 
Vi _- Fan cae Spi = é 
3 H My 


30 


Ais MMe in which the mass, Mr,is the mass of the target atom. 


Mi +Mr 

Define ET = yy (4 et ee) sp (7) 
Define V= VER) 

ET (8) 
Then &, is the root of . 

|\-~-%+-S=0 (9) 
The root of this equation was found using the Newton-Raphson method 
which is ik; = Rk, = £ CA») (10) 
#"(Rn) 

where = _ Bee ics 11 

f (R) R = (11) 
and ‘ton -  VESIGDY . (12) 

£ (R) - t & 


in which DV is defined by 


nve Me Cres apes cre”) (13) 
then dk ET 


_ _ (RR -veR, - © 
Raer = Re V = RyDV +207R, pe 


The integral, U, was evaluated using the Gauss-Mahler Quadrature 


(1) 


formula as employed by Bates, Cook and Smith The result is 


U= oh, EF, (cosse)acost +(CosH)e(Cosh) +(cosfh)g(cos #) (15) 
2,2 - '/2 
gen.) =[1 - Yi - Sx] a6) 


where 


Q 


in which the 
Xi = Cos# , cosh and cossf (17) 


To calculate |22 differentiate equation 2 witth respect to@ 


Fa = j2egu + 2U | (18) 
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where 


3h =F osegge (cos) RMD os gigeeer 


and the X. are as defined previously. 


dR 
ar is obtained by differentiating 


3 
the result is 


i oe 
d@ 9-2 O*/p, -R,DV+V 


aU = (costes + (cos BRR) (oss gate ©) ede 


in which 


2 

-* 
9 9G MX = A - VROX; = @* x? DV. — 
OR, 2 Ro R 


where the X, are as defined previously, 
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EXPLANATION OF PROGRAM 


STATEMENT INDEX 


1-5 Places and prints column headings for print out. 
6-8 Dimension statements. Self explanatory. 

9-23 Do loop for graph title. 
24-29 Constants for the potential used. See equation 


No. 6, Appendix 1. May take on any values in 
accordance with potential used. 


30 Sets the initial value for the number of 
iterations. See statement index 74 and 75. 


31 Sets the initial value of the impact parameter 
such that RHO(1) = 9. See statement index 60 


32 Sets the initial value of theta, the scattering 
angle in C. M. Coordinates. 


33 Sets the initial value of the energy of the 
hydrogen molecular ion in the LAB. Coordinate 
system, EV, such that EV(1) = 1000ev. See 
statement index 37. 


34, 39, 40, 41, Give the draw program the proper parameter to 
42, 43, 44, 45 draw 10 curves on one graph. 
35 Begins outer do loop for 10 different energies 


defined by statement index 33 and 37. 
36 Subscript relationship to achieve desired energies. 


oy Energy defining relationship. May be varied for 
desired energies. 


38 Sets the initial value of R for the iteration 
to find the root of equation 4, Appendix 2. May 
be based on b, the collision diameter. (Pagel 9) 


47-57 Label program for the labeling of curves with 
appropriate energy. 


58 Begins inner do loop for 101 values of the impact 
parameter, RHO, as defined by statement index 31 


and 60. RHO is in units of ae 


a9 Subscript relationship to achieve desired values 
of RHO. 
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60 


61-63 


64 
65 


66 


67 


68 


69 


70 


71 


72-73 


74 


ri, 


76 


77 


78 


79-61 


82 


83-85 


86-100 


Defining relationship for RHO. May be varied 
as desired. 


Evaluates the exponential terms of the potential 
used. 


Mass of the projectile particle in a. m. u. 
Mass of the target particle in a. m. u. 


Evaluates the relative energy factor of equation 
4, Appendix 1 in terms of EV. 


Evaluates V as defined by equation 8, Appendix l. 


Evaluates equation 11, Appendix 1, for the 
iteration. 


Evaluates DV as defined by equation 13, Appendix 
1, for the iteration. 


Evaluates equation 12, Appendix l. 

Computes a new value of R for the iteration. 
Compares the last two values of R. 

Counts the number of iterations. 


Stops the program if finding R requires over 
10000 iterations. 


Replaces R with the last value computed by 
statement index 71. 


Repeats the iteration if statement index 73 has 
not been satisfied. 


Continues the program if statement index 73 has 
been satisfied. 


Insures that the proper value of R has been used 
to compute V and DV. 


Stops the program if R is less than RHO and exits 
the inner do loop if R = RHO. 


Evaluate the cosine functions of equation 15, 
Appendix 1. 


Evaluate the g (X, ) as defined by equation 16, 
Appendix l. 
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102 


PG3-105 
106 
107 
108-113 
114 
115 
116 
Py 
118 
Try 
120 


ae 
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124-127 


128-132 


Evaluates the integral U as defined by equation 
15, Appendix l. 


Computes theta the scattering angle in the C, M. 
coordinate system for the given energy and impact 


parameter by equation 2, Appendix l. 


Computes the O9%)4e as defined by equation 21, 
Appendix l. 


Evaluates 9U/9@ as defined by equation 20, 
Appendix l. 


Computes dR./ye@ as defined by equation 22, 
Appendix l. 


Computes the BIKY, 2 as defined by equations 24, 
Appendix l. 


Computes aV/5R, as defined by equation 23 
Appendix l. 


Computes dU/do as defined by equation 19, Appendix 
Ibe 


Computes }d®/0} as defined by equation 18, Appendix 
1 


Computes the differential cross section for 
scattering of protons at the previously computed 
angle by equation 1, Appendix l. 


Computes A pas given on page 16. 


Prevents dividing by zero in statement index 122 
when theta is 180. 


A dummy statement to give a value of A @® for 
theta = 180. 


Skips statement index 122 when theta = 180. 
Evalue ae given on page 16, 

] 
Evaluates sumce given on page 16, 

2 


Compute the percent difference in velocities 
after scattering as described on page 24, 


Convert the scattering angle, theta, to LAB 
coordinates and degrees. 
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133 


134 


135-136 


ES 7) 


138 


139 


140 


141-142 


143 


144-150 


Gives the inter atomic spacing &. 


Computes the value of the orientation probability 
as defined on page 


Limits the value of P to l. 


Computes the differential cross section for 
H, scattering as given on page 13, 


Gives the proton energy for printing when that 
is of interest. 


Prevents the log. function from having a negative 
or zero argument in statement index 140. 


Computes the logrithum of the crosssection for 
the draw program. 


Limits the values of the scattering angle printed 
to those between 2° and 90°. 


Prints results for each value of RHO assigned. 


Self explanatory. 


41 


INITIAL DISTRIBUTION LIST 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library 
Naval Postgraduate School, Monterey, California 


Defense Atomic Support Agency, Dept. of 
Defense, Washington, D. C. 20301 


Professor Robert L. Armstead 
Department of Physics 
Naval Postgraduate School, Monterey, California 


LT Eugene Paul FEIST, USN 
Engineering Department 
Naval Academy 

Annapolis, Maryland 21401 


LT Tom Bush 
SMC # 2107 
Naval Postgraduate School, Monterey, California 


Dra. vs cCook 
Stanford Research Institute 
Menlo Park, California 94025 


Professor 0. Heniz 
Naval Postgraduate School, Monterey, California 


42 


No. 


Copies 


20 


ke 


Securit Classification 





DOCUMENT CONTROL DATA - R&D 


(Security claesificetion of titie, body of abstract and indexing annotation muaet be entered when the overall report ie claeeified) 


1. ORIGINATING ACTIVITY (Corporate author) 
Naval Postgraduate School UNCLASSEFIED 






25. GROUP 


Monterey, California 93940 

3. REPORT TITLE 
Classical Calculation of the Differential Cross Section for 
Scattering of Hydrogen Molecu@ar Ions From Heavy Atoms 


4. DESCRIPTIVE NOTES (Type of report and inclusive dates) 


Ga. CONTRACT OR GRANT NO, 9a. ORIGINATOR'S REPORT NUMBER(S) 


é 


&. PROJECT NO. 


9b. OTHER ek NO(S) (Any other numbers that may be aselgned 
thie repo 


He rn oper ae OP eee, 


11. SUPPL EMENTARY NOTES 12. SPONSORING MILITARY ACTIVITY 





13. ABSTRACT 


The primary problem in the calculation of the differential 
scattering cross section for hydrogen molecular ions is the derivation 
of an approximate expression for the probability that the hydrogen 
molecular ion will be oriented, prior to impact, in such a way as to 
avoid disruption during scattering. Using a simple model, an 
approximate expression for the orientation probability is obtained. 
Differential cross sections for the scattering of hydrogen molecular 

: ions from carbon and krypton are calculated using classical 
scattering theory and the derived orientation probability. Impact 

! energies of 100 to 1000 electron volts are considered. In addition 
differential cross sections for protons scattered from carbon and 
krypton are presented. 


DD ."°%.1473 __INGLAS ST rrEp 


43 Security Classification 


a! 
4 


~ Security Classification 


14 
KEY WORDS 


Scattering Cross Section 





Hydrogen Molecular Ion 


Orientation Probability 


FORM | , 
DD 1 NOV 31473 (BACK) UNCLASSTFIED 


S/N 0101-807-6821 aa Security Classification A-31409 

















thesF2526 a 
Classical calculation of t 


EY KNO 


DU 


WME AL 


Pinte 


™  DUDUULEY KNUA LIBRARY 


F. 
2%, 
! 





